Strongly Inelastic Granular Gases 
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The expansion of the velocity distribution function for the homogeneous cooling state (HCS) in a 
Sonine polynomial series around a Maxwellian is shown to be divergent, though Borel resummable. 
A convergent expansion for the HCS has been devised and employed to obtain the HCS velocity 
distribution function and (using it) the linear transport coefficients for a three dimensional monodis- 
perse granular gas of smooth inelastic spheres, for all physical values of the coefficient of normal 
restitution. The results are in very good agreement with findings of DSMC simulations. 

PACS numbers: 05.20.Dd, 45.70.-n, 47.45Ab, 45.70.Mg 



The success of granular hydrodynamics in describing 
granular gases is by now undeniable [l|. However most 
favorable agreements between theory and experiment (or 
simulations) have been restricted to mild degrees of in- 
elasticity, i.e. coefficients of normal restitution, a, that 
are larger than about 0.7. Furthermore, DSMC simula- 
tions seem to agree with theoretically computed values 
of the transport coefficients [2fl for 0.7 < a < 1 but def- 
initely not for a < 0.6. This is understandable for the 
comparisons which provides accurate results for 

the near-elastic case, but less so for [|[ which is formally 
correct for all values of a. The reasons for this situation 
and its resolution are explained below. 

Consider a dilute gas of inelastic smooth and homo- 
geneous hard spheres of mass, m, diameter, a, number 
density, n, and coefficient of normal restitution, a. The 
Boltzmann equation for this case is [H, 0, Q : 



df(r,v u t) 

at 



vrV/(r, Vl ,t)=B (/,/), 



(1) 



where / denotes the velocity distribution function, Vj 
denotes the velocity of a sphere, and r the position of its 
center of mass. The collision term B(f, /), is given by: 



1 



~ 0/1/2 _ /1/2 



B(f,f) = a 2 f dv 2 dk(k-v 12 ) 

Jk vi2 > 

(2) 

where fa = / (r, Vi,t), fl = f (r, v£, t), k is a unit vector 
along the line of centers of the colliding spheres at contact 
and primes denote precollisional velocities. For any set 



of vectors, {A,}: Aj 



Aj. The binary collision 



between spheres "i" and "j" is given by: 



1 



(kVy)k. 



(3) 



Two systematic perturbative (generalized Chapman- 
Enskog) expansions for solving Eq. (Q]) have been devel- 
oped. Ref. [|| proposes a double expansion in the gra- 
dients (or the Knudsen number, K) and the degree of 
inelasticity e = 1 — a 2 about a local Maxwellian, Jm- 
Ref. [1] employs a gradient expansion around the local 



HCS distribution function, / W . Following the Chapman- 
Enskog method, it is assumed in both approaches that / 
depends on time through the hydrodynamic fields: the 
number density, n (r, t) = J f (r, v, t) dv, the (granular) 
temperature, T (r, t) = / / (r, v, t) w 2 dv = fv 2 (r, t), 
and the velocity field, u (r, t) — J f (r, v, t) vdv. Define 
the peculiar velocity by: V = v — u. It is convenient 
to nondimensionalize velocities by the thermal speed vo', 
in particular define: c = V/t>o, and the dimcnsionlcss 

distribution function, / (r, c, t) = —f (r, v, t). 

Consider first the HCS. In the absence of gradients, 
Eq. (JlJ admits a solution, /^(c), where c = |c| (this 
'scaling solution' is known to be a 'long-time' limit of the 
homogeneous solution for a non-HCS initial condition, 
but this detail is not relevant here). It follows that for 
the HCS: 



/<°> (c)dc=l 



c 2 / (0) ( C )dc = ^. (4) 



Eq. |T]) becomes 



where 



eW 1 



W 



ci d 
YdcT 



f 0) = B(r>,f 



f(0) f(0) 



dddcacJa/W ( Cl )/W ( C2 ) , 



(0) 



(5) 



(6) 



andfl (/ (0) ,/ (0) J is the non-dimcnsionalizcd Boltzmann 

operator ((21). It is common to represent /^°^ in terms 
of a Sonine polynomial expansion 0, / < - - ) (c) = 
/jVf(c)$(c), where /m (c) = ir~%e~ c is a Maxwellian 
distribution, and $ (c) = J2^Lo a p^l {° 2 ) ■ The substitu- 
tion of this expansion, truncated at p = 2, in the Boltz- 
mann equation yields values for the coefficients up to a-i-, 
see HQ. Ref. presents indications that this expan- 
sion diverges and the lack of convergence is attributed to 
the well known 0] exponential tail of the HCS (at large 
speeds): £^ ^ (c) ~ e~^ c = f^ u (up to an algebraic 
1). 



prefactor [9|). Indeed, a direct calculation of the expan- 
sion of $, for /(°) replaced by 7fi f^ u (namely the contri- 
bution of the tail to the expansion), i.e. an expansion of 



2 



y-lds |ld: a p = % £^ ^1 (^) 

and for large p, a p ~ {-f p^ (^f P+Z e"W, 
i.e. the series is divergent, though asymptotic and 
Borel resummable. Furthermore, using the method ex- 
plained below we have computed the ratio — ;z7~jrn T' 
vs. p for the correct distribution of the HCS, where 
z = ML, for a = 0.1 and z = 0.544 and find (us- 
ing N s — 40 polynomials, see Fig. [I]) that it is es- 
sentially constant for p > 5. It is also interesting to 
consider the sum (in which the asymptotic form of a p 
is substituted for a p ): S = J2pto(-) P z 2p P ] -S{(c 2 ). 



t-Wq+T) /2eW\ 2 <l+ 3 




FIG. 1: — tt-ttt^ 1 vs p. a = 0.1, z = 0.544 and N s = 40. 
The descent beyond p = 10 is a result of the truncation error. 



polynomials is divergent, though Borel resummable. 

Rather than using Borel resummation to obtain the 
distribution function of the HCS (and later the linear 
transport coefficients) we propose a method that pro- 
duces convergent series for the distribution functions. To 
this end define the following expansion for the HCS: 



/ (0) (c) 



^-fe-^ 2 ^ a (7)^|( c 2 

p=0 



(8) 



where 7 > is a constant, i.e. we expand the HCS around 
a Maxwellian corresponding to a 'wrong' temperature. It 



can be shown that the coefficients a p ^ of e c * ' 



tan satisf y 



,(7) 



x 



— " — 3" (1^7) ' nence a convergent expansion 



is expected (and obtained) for 7 < 1/2 (in practice we 
chose 7 = 0.4). Upon substituting Eq. ([5) into Eqs. (0]- 

IHl), multiplying Eq. ([5]) by S± (c 2 ) and integrating over 

2 

the velocity one obtains: 



p=0 



,(7) 



(7-l)\t- 



p,q=0 



p=0 



2' 



(9) 



B Npq aWaW =eWj2 M Np a$\ (10) 



This sum can be Borel resummed (the Borel sum is 
denoted by Sb) upon replacing pi by J Q e~*i p di and 
exchanging the orders of summation and integration: 

Sb = J™e- t J2 P =o(-tz 2 ) Ps i ( c *) d< - Next - the sum 
under the integral can be calculated by inspection noting 
that the generating function of the Sonine polynomials is 

Er= ^w = ( i - s )" m " le "^ : 



Sr. 



dt. 



(7) 



where 



M 



N P 



Qr, 



p,q=0 



2 / 7 _iyr(p + §) 



y/TTJ2 \ 7 



(11) 



Using the Laplace method one obtains from Eq. (J7J) a 

result which is practically identical to e° 2 f^u- Sb ~ W p , 

r ( 371 c eW ) 

iie %/3tt C ; where C is a constant. Thus we have 

c 

shown that the standard expansion of the HCS in Sonine and 



871-i 



d Cl dc 2 e^( c M C ? 2 S? (c?) S« (c 2 ) 



SjVpg = 7T 3 



dcidc 2 dk (k • Ci 2 ) S 



JV 



k ci2 > 



i_ e -7(c' 1 2 +4 2 ) 5 P ^,2) 5 | ( c /2) _ e -7(c?+c 2 ) 5 P ^2) ^ 



(12) 



When the above sums are truncated at iVg, one ob- tains a quadratic algebraic system in Ns + 2 variables: 



3 



W, ag 7 ^ , . . . , ajj* . The pertinent prefactors are obtained 



as follows. Define the "super-generating" function [ill ] : 



I (h,b 2 ,b3,b4,,x,y,z) = / dcidc 2 dk (k • c 12 ) e F , 
where 

F = bicf + b 2 c 2 2 + b 3 c? + he? + | ( Cl + O 2 

+ |(c 1 +c' 2 ) 2 + ^(c 1 +c 2 ) 2 . (13) 

This integral can be carried out to yield: 



I {b 1 ,b 2 ,b 3 ,b 4: ,x,y,z) 



AiAJ (Ai + A 3 ) 



(14) 



where Ai 



_ A 2 y+z (b 1 -b 2 +b 3 -b 4 +2x) 2 



4 2 4A 2 



2(x + y + z), X 3 = ^(b 3 + b 4 +^ 1+a 

Next, using Eqs. (fTZl [T4")) define a generating function 
from which B^ pq can be obtained by taking derivatives 
(a symbolic processor is used for this purpose): 

00 

B(r,s,t)= r N sHW Npq = 

N,p,q=0 

7T -3 (1 - r)~^ (1 - s)~5 (1 _ t)-§ x 



1 



r s t 
,0,- : + 7 



2 v 1 _ r ' ' 1 — s 
r s 



1 -t 



7,0,0,0 



1 — r 1 — s 



+ 7, 



1 - t 



7,0,0,0,0,0, 



The coefficients W pq , M^ p and Q p are obtained in a simi- 
lar fashion. The solution of the Ns+2 algebraic equations 
for {a n } now enables one to obtain f(°>. Fig. [2] depicts 
the convergence towards the exponential tail for a = 0.1 
through a plot of — ^ ln/(°) vs. c (which should be con- 
stant for an exponential function) for N$ — 10,20,40. 
The predicted prefactor in the exponential is and 



using Eq. 



-JjL m /(o) (3) _ 3.686. 
The coefficient a p (= 



it equals 3.679, while, e.g., for c = 3, 



) is the projection of / (°) on 



S? 



/ dcS{ (c 2 ) /(°) (c). Fig. presents 



' ~ 2r(n+f) 

a comparison of the correct value of a 2 with the result 
obtained from a truncation at p = 2 



«2 



16(1 -a) (l-2a 2 ) 
81 -17a + 30(1 -a) a 2 



(15) 



Next, we apply our results to the theory of Ref. 
The linear order (in gradients) in the Chapman-Enskog 
expansion for Eq. © is /W = ^(V)-VlnT + 
B ( V) -Vlnn + Cij (V)diUj, where an overlined tensor 




5 2 2.5 3 3.5 



FIG. 2: -4Mn/ w vs. c for a = 0.1. The Sonine sum was 

truncated at N s = 10 (dashes), N s = 20 (dots) and Ns = 40 
(solid line), to show the convergence towards the tail. 
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FIG. 3: A plot of Eq. (ll5|l (solid line) compared to the con- 
verged value of 02 (asterisks) vs. a. 



denotes its symmetric and traceless part. The zeroth or- 
der cooling rate satisfies the following equalities 0]: 
C (0) = -c\ (0) InT = f n<r 2 v W where W is the solution 
of Eq. ((6]). Define the linear operator 

Lh = -B(f^,h\-B(h,f^y (16) 

The functions A(V), B (V) and (V) satisfy Q: 



c + £pj B = c (0 U-v/(°)-|a v / (0) , 



((0) 

~2 

(•(0) 



(17) 



where £ = L- C [0) Td T - ^- and the notation of ref. 
is employed. 

At this order, the pressure tensor, P.^ = mj dvViVjf, 
and the heat flux, q = m/2 J dvV 2 V/, are: 



P-P = md k ui I dvCkiViVj = -2r]diUj, 



4 



<zf 3 = j J &W 2 V t {A k d k \nT + B k d k \nn) 

= —k T dilnT ~ fi n dihin, (18) 

thereby defining the shear viscosity, 77, the thermal con- 
ductivity, k, and the diffusive heat conductivity, fi. De- 
fine [i| the 'reduced' coefficients 77*, k* and /i* as the 
transport coefficients normalized by their respective elas- 
tic values evaluated at the lowest order in the Soninc 
expansion: 77* 

m _ 5mv 



^ k* = and n* = Tt ^ where 



and kq 



75k B vo 



We have solved 



Eqs. (|17[) using a Sonine polynomial expansion with N$ 
up to 10 (unlike the tail of the distribution, the trans- 
port coefficients depend on low order moments of Z^ -*), 
for which convergence was attained. Figs. |U [5] and [5] are 
plots of 77*, k* and /i*, respectively as a function of a. 
The agreement between DSMC simulations ('o') and the 
computed coefficients ('*') is excellent, except in the case 
of 77* . For comparison we have plotted the analytic re- 
sults obtained in the first Sonine approximation 
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FIG. 4: The reduced shear viscosity, 77*, as a function of a. 
Asterisks denote the theoretical values and circles represent 
the results of DSMC simulations. The solid line is the ana- 
lytical result obtained in the first Sonine approximation [j. 
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FIG. 5: The same as in Fig. @ for k*. 
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FIG. 6: The same as in Fig. (0J for /j*. 



In summary, we have shown how a convergent expan- 
sion can be successfully employed to obtain transport 
coefficients for all physical values of the coefficient of 
restitution. The direct physical relevance of these re- 
sults depends, among other things, on the importance 
of precollisional correlations for strongly dissipative sys- 
tems (which may require going beyond the Boltzmann 
equation), as well as that of collapse events. 
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